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Conditions for smooth cosmological models are set out and applied to inhomogeneous spherically 
symmetric models constructed by matching together different Lemaitre-Tolman-Bondi solutions to 
the Einstein field equations. As an illustration the methods are applied to a collapsing dust sphere 
in a curved background. This describes a region which expands and then collapses to form a black 
hole in an Einstein de Sitter background. We show that in all such models if there is no vacuum 
region then the singularity must go on accreting matter for an infinite LTB time. 

- - - I. INTRODUCTION 

O 

Recently we and others §M3 have had cause to use matched spherically symmetric solutions of the field equations 
CN ■ to model fractal || distributions of matter in cosmological situations. Single spherically symmetric models (see [Q for 
a recent example) and matched ones have long been used in cosmology (see [pi for an extensive and comprehensive 
review). They are remarkably rich in structure but there are many subtleties in their application to cosmology [|6|-gJ. 
The purpose of this paper is to clarify and unify existing results on the construction of smooth models. The topic has 
important implications for the modelling of structures like voids or the formation of black holes in curved backgrounds. 
The new results include matching between exact solutions in the Kantowski-Sachs family and the Lemaitre-Tolman- 
Bondi (LTB) solutions and to the occurrence of centres. It is not always realised that a spherically symmetric model 
need not possess a centre |1 . Here an extended list of the possible types of centre is given. Also all possible regular 
composite dust models are given and they are classified into four classes according to the number of centres they 
allow. One of the classes in which the spatial sections have the topology of a 3-torus appears to be new. 

This work complements that by Lake |10| and by Fayos et al. fLlfl . The work by Lake is more theoretical. That 

by Fayos, Senovilla and Torres is more general and very interesting geometrically but we find it less easy to see some 

subtle problems such as the one we illustrate in the final section. It has long been accepted that a black hole in 

an Einstein-de Sitter model can be constructed using a Schwarzchild interface. What had not been shown but is 

^D ■ demonstrated here is how a black hole evolves from an expanding region in an Einstein de Sitter background. The 

model is a modification of one produced by Papapetrou |20| which is non-physical because the singularity takes an 

QT 1 infinite LTB time to be completed. In section 7 it is shown that this infinite time is an inevitable result of the matching 

J conditions if no vacuum region is included, a result which appears to be not widely known. 
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II. SPHERICAL DUST SOLUTIONS 

In comoving coordinates x a — {t, r, 9, <f\ the spherically symmetric metric can be written [|jj 

ds 2 = -dt 2 + X 2 (r,t)dr 2 + R 2 {r,t)dn 2 , (1) 

and the dust 4- velocity as u a — Sf. Here dfl 2 = dO 2 + sin 2 9 dej) 2 and R > with R = only at a centre. The 
coordinate ranges are; r > 0, < 9 < ir, < cf> < 2tt and t > T(r) (see the solutions below for a definition of T). The 
dust energy momentum tensor is T fc a = fiu a Ub, where /j, is the proper matter energy density. 
In order for the Einstein tensor to be well defined we require that 

X ^ ± R; R is C 2 in r and t; X is C 2 in t and C 1 in r. 

We call such behaviour regular. Regions of regularity may be joined together to form composite space-times in which 
the differentiability conditions hold only piece-wise. More details are given below when we discuss the matching. To 
set the notation and for later use we first write down the first integrals and write out the solutions |l2,0,m|. We will 
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call the solutions LTB models and refer to the coordinates as LTB coordinates. We use units such that G = c = 1 
and the notation, overdot for u a d a = d/dt and prime for d/dr. If R' ^ 0, then the Einstein field equations reduce to 

x = ^Wr (2) 

where E(r) is an arbitrary function of integration, and 

te-^ + E, ( 3 ) 

with M(r) a second function of integration. The proper density is given by 



M ' , A \ 

/i= 4^^' (4) 

The equations (@) - (0) have five solutions: 

[si] for {E = M = 0}, R=-T(r), 

[s2] for {E > 0, M = 0}, R = VE(et - T), 

[s3] for {E = 0, M > 0}, R = {9M/2) x / 3 (et - T) 2 / 3 , 

[s4] for {E > 0, M > 0}, R = f-(cosh7 ? - 1), 

sinh?] - ?/ = (et - T)E 3 ^ 2 M- 1 1 
< r\ < co, 

[s5] for {E < 0,M > 0}, i? = M(cos?7 - 1)-ET\ 

77 - sin?; = (et - T)\E\ 3 / 2 M-\ 
< j] < 2n. 

Here, T(r) is a further function of integration and e = ±1. The solutions si and s2 are locally Minkowskian. The 
surfaces {et — T = 0} are spacelike and singular. Following Bondi 0, M is the relativistic generalisation of Newtonian 
mass and ^E is the total energy. 

The case R' = leads to an inhomogeneous generalisation of the Kantowski-Sachs metric llq] Gq] which does not 
have centres of symmetry in the hypersurfaces t = constant. It is given in our notation by 



[s6] 



J R = M(l-cos7y), (5) 

r)-smr] = M- 1 (et-T), (6) 

77 sin 77 



X = A{r) sm7? +B( r ) 
v 7 1-cost7 v 7 



1- 



(7) 



2(1 — cos 77) 
where A and B are integration functions, Mis a constant, and < 77 < 2ir. We label this solution s6. 

III. JUNCTION CONDITIONS 

Our aim is to establish a class of cosmological models which do not have shell-crossing or surface density layers, 
but which can be constructed by matching together different LTB solutions from the set si to s6, so we impose the 
Darmois conditions Km on the junctions. For a comoving space-like junction r = constant and dust, these require 
that the first and second fundamental forms are continuous across the junction (interface). The unit normal is given 
by n a — l^l^o) an d the first fundamental form is the metric intrinsic to the interface, i.e., 



h a b = gab - n a n b 

= diag(-l,0, R 2 , R 2 sin 2 9), 

and the second fundamental form is the extrinsic curvature which is given by 

Kab = h c a htV d n c 

,. / „ RR' RR' . 2 
= diag 0, 0, -r^rr, tttt sir ■ 



|X| ' \X\ 

From these it follows that the necessary and sufficient conditions for matching are that 

R is continuous in r, 



(8) 
(9) 



(10) 
(11) 

(12) 



Rf_ 

\X\ 



is continuous in r. 



(13) 



The nature of the problem changes in the non-comoving case 1 16 JL7| . The conditions at the junction must be satisfied 
through some range of values of r. Explicitly, if the spacelike non-comoving boundary surface is given by 



r — g(t) — constant, 



where for convenience we choose P such that 



P 



dg_ _ 

dt \X\VP 2 -l' 



then the unit normal to the surface is given by 



n a = P, 



= ,0,0 



\x\VP 2 -i' 

and n a n a = — 1. The unit tangent to the surface with 6 — <p — constant is 

P 



m a =(^2-1 o,0 



The intrinsic curvature of the surface is given by 

P 2 - 1 -P\X\y/P 2 - 1 







-p\x\Vp 2 - 1 p 2 x 2 o o 





R 2 i? 2 sin 2 ( 



and the extrinsic curvature is 



K ab = 



( (P 2 - l)Fi 



-F^XIPVP 2 -! 

-Fi\X\P^/P 2 -1 P 2 X 2 Fi 

F 2 











F 2 suv z ej 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



where F\ = ^; + P^ an d Fi = PR^rr ~ m /t-i i an< i a U f ne ordinary derivatives are taken along the paths {9, </>} = 
constant in the tangent space to the hypersurface. On the hypersurface, one coordinate is surplus because of (|l4J). 

We will not deal with non-comoving junctions further here except to mention two points. First that at the junction 
at least one of E, M and T must be constant for a range values of r. This follows from the fact that they have to be 
continuous by the Darmois conditions and at least one of them has at most one value in common (or as a limit for 
s6) for allowed[J matches between the solutions si to s6. Thus in the non-comoving case there will be an interface 



1 see later 



region. Second, that Krasinski |L8| suggests that non-comoving boundaries could find an application in the formation 
of structure. The idea is that it may be possible to use them to allow incoming matter to augment the condensations 
discussed by Bonnor |1S[] , 

From here on we restrict attention to junctions where there are no surface layers or interface region and so they are 
necessarily comoving and the metrics on both sides are of the junction are determined by the integration functions 
E(r), M(r) and T{r) (and A(r) and B(r) in regions where E{r) = —1). For the spherically symmetric dust models 
the junction conditions (O) imply that 

M(r)(> 0), E{r)(> -1) and T(r) are continuous. (20) 



IV. REGULARITY REQUIREMENTS 

In this section the physical requirements to be imposed on the metrics are made more explicit and justified. It 
follows from the conditions on the metric and the matching conditions across spacelike surfaces that the metrics si 
to s6 are at least C 2 in t. Thus we are only concerned with the behaviour with respect to the radial coordinate r. We 
start by requiring that lim±A ^ everywhere including across the junctions. This condition is important because it 
enables us to express the continuity properties of physical quantities unambiguously through their differentiability in 
r, and prevents shell crossing. To ensure that there is not a curvature singularity as R — > 0, we require that M/R 3 is 
finite everywhere except trivially at the spacelike singularity r — > 0. Finally we require that at a centre the shear 



^=3 



X R 
X~R 



x diag(0,2,-l,-l) 



go to zero to maintain spherical symmetry. If it is not zero, then the eigendirections of the shear tensor will violate 
spherical symmetry at the centre. 

In summary, we require the following properties to hold in our matched spacetimes: 

Rl The junction conditions (|12,0) hold. 

R2 lim ± A ^ 0. 

R3 The shear tends to zero whenever R — > 0, i.e. at a centre. 

R4 M/R 3 remains finite as R — > (except trivially at the spacelike singularity r — > 0) . 

R5 The metrics si to s6 are regular within their domains, i.e. between junctions. 

We will now describe the implications of these conditions within the domains of metrics. It follows from the field 
equation (0) that 

R' can change sign only at values of r for which E(r) = — 1, (21) 

i.e., in s5. For si to s4, R' > or R' < throughout. For metrics si to s4, differentiation of the exact solutions for 
R, and examination of the asymptotic behaviour for large and small t, yield 

± R' > => {±M' > 0, ±E' > 0, ±T' < 0}. (22) 

For s5 in the region in which E(r) =/= — 1, i.e. where R' does not change sign, differentiation of the exact solution for 
R with respect to r and investigation of the result as r\ — > 2-7T, and assuming [i > 0, gives 

±i?'>0^ |±M'>0, ( E'M - -M'E + — 1— ! — J > 0, ±T'<o|. (23) 

Relations (g2|) and (|2^) are the Hellaby and Lake || no-shell-crossing conditions. 

For values of r at which E(r) = — 1, the condition 1/0 and finite implies that we must have 

R! ^ 



and finite as r — * r* where E(r*) = — 1. Again, differentiation of the solution for R with respect to r gives 



lim r ^ r » 



M' 



lim,.^ r . -f, „ lim r ^ r . L 



n+E 



l + E 



must be finite, and at least one must be non-zero. For a non-zero density at r* 



M' 
/T+E 



must be non-zero. 



More generally, // = in si and s2 from their definition. In s3 and s4, fi is only zero if M' = 0; otherwise \x is 
finite for all r since B! ' jM' ^ 0. For s5 it follows from above that R' JM' ^ 0, and hence the density is finite even 
where E{r) = —1. For s6 the density vanishes if B — and A^O, otherwise it is finite and positive if and only if 

,• B 1 

hm — > — . 

± A ~ 7T 

We now consider the conditions for a centre (R = 0) to exist or be attached to a solution. Only comoving centres are 
possible and they may only join to solutions si to s5. A list of the possibilities is given in table 1 |16|. 



Soln. 


Behaviour of E, M and T 


Kinematics 


Example 


(si) 


(i) T - 

(ii) limT' finite, nonzero 


e = o 


T = -r 


(s2) 


(i) E - 

(ii) lim(_E 1/2 E') finite, nonzero 

(iii) ET'/E' -► 




£ = r 2 
T = 


(s3) 


(i) M -> 

(ii) lim(Af ~ 2/3 Af) finite, nonzero 

(iii) MT'/M' -> 


e -» 26T- 1 

47r^i — » |r" 2 


Af = r 3 
T = 


(s4) 


(i) E 3/2 /M -» 0, M -» 

(ii) lim(Af ~ 2/3 Af') finite, nonzero 

(iii) lim(Afr'/M') = lim(Af 1/3 £'/Af) = 


9 -> 2er" 1 
47r/i — » |r" 2 


£ = r 3 
Af = r 3 
T = 


(i) E a/2 /M -> +oo, £■ -> 

(ii) lim(_E -1 ' _B') finite, nonzero 

(iii) lim(ST'/S') 

= lim [iT^M'JS'- 1 log(E 3 / 2 M- 1 )] = 


4-Kfi — » 


£ = r 2 
M = r 4 
T = 


(i) E 3/2 /M -* a > 0, M ->■ 

(ii) lim(_E~ 1 M') finite, nonzero 

(iii) E- X ME'/M' -> 2/3 and MT'/M' -» 


6 — * 3ea sinh rj/(cosh 77 — 1) 
Anfj, -» 3a 2 /(cosh 77 - l) 3 


£ = r 2 
Af = r 3 
T = 


(s5) 


(i) |£| 3/2 /Af -> 0, Af -» 

(ii) lim(Af ~ 2/3 Af') finite, nonzero 

(iii) lim(AfT'/M') = lim(Af 1/3 £'/Af') = 


-> 2CT- 1 
A-Kfj, — > §r~ 2 


E = -r 3 
M = r 3 
T = 


(i) |£| 3/2 /Af -► a > 0, Af -» 
(ii) lim(_E _1 M') finite, nonzero 
(iii) E~ X ME'/M' -> 2/3 and MT'/M' -> 


6 — » 3ea sin 77/(1 — cos?;) 2 
An\\i — » 3a /(l — cost;) 


£ = -r 2 
Af = r 3 
T = 



TABLE I. Central Behaviour 



In the table, results labelled (i) arise from the behaviour as R — > 0; those labelled (ii) derive from the requirement 
that X ^ 0, and (iii) are a result of the shear vanishing. For each case, M/R 3 — > |7T/x, which is the Newtonian limit. 
Simple examples for which the centre lies at r = are given in each case for illustration. The expansion rate at the 
centre, 



e 



'R 



X 
X' 



is listed for each solution. The central behaviours listed for s4 and s5 j 
the calculations involved to derive the results in the table is given in pi 



;eneralise previous results. An illustration of 



V. MATCHING OF SOLUTIONS 

In this section we reach the core of the paper. We examine matching across comoving space-like surfaces between 
solutions si to s6 to form composite models. The sign of R' cannot change across these interfaces since E ^ — 1 on 
them - except on interfaces between s5 and s6. However note that R' = in s6. Solutions si may not be matched to 
any others because they are cosmological and it has R = 0. Solution s2 does not match to s5 since (E ~ * 0,M — > 0) 
forces R — * in s5. From the properties that characterise s6 as an LTB model (see section 5.5), it follows that it only 
matches to s5. There remain just five physical types of junction. 

A. (a) Matching s2 to s4 

The s2 (interior) side of the junction is unconstrained by the matching. Approaching the junction from s4, M — > 0, 
E > and M' > (M' is continuous in s4) since M > in s4 and M = in s2. Hence R' > and R is increasing 
in the direction s2 to s4. Denote the value of r at the junction by r* . Then as r — > r* , r/ — > oo because M — > and 
E > in the exact solutions (s2) and (s4). It follows that 



X 



VTTe 



M' [E 3 / 2 \ , 1/2 E't 
-E l0 4-W)- TE/ + 2EV- 2 



and hence that X satisfies R2 if lim( s4 ) T',lim( s4 ) E' and lim^) M'lnM are finite and at least one is non-zero. 
On the (s4) side, the density reduces to 



47T/i 



r 2 log 



E 2 / 3 T't 2 \ E 3 / 2 T't 2 E^^E't 3 



M 



M' 



2Af 



(24) 



which goes to zero asr^r*. Note that \i = in (s2). 

Since on the (s4) side M — > as r — » r* and M'M is finite, we must have M' — > and hence 



X 

X 







On the (s2) side 



1ET' 



X 



2M' 



E 1 / 2 E> log M 



if E' -» 0, 

otherwise. 



A r - 2ET'/E> ' 
and on both sides -R/i? — > e/r. From this it follows that the shear remains finite on both sides of the junction. 

B. (b) Matching s3 to s4 



The s3 side is unconstrained by the matching. Approaching the junction from within the s4 region, E — > and 
A/ > 0. Since E > in s4, we must have E' > in a neighbourhood of the junction in the s4 region. Also, from the 
explicit solution for s4 we have 

sinh?7 - 77 = tE 3/2 M- 1 -» 0, 



because E — > and M > 0. Thus ry sa (6r/M) 1 / 3 £ ,1//2 near the junction on the s4 side and E' > 0, since E > in 
the (s4) region. Hence R must increase in the direction (s3) to (s4). On the (s4) side, 



X->M 



1/3 



1/3 



,( r 2 Y /J r (Uiy> E' f(6rY 



T \^) 



40 



M 



1/3 



V6M 2 
and so X is finite and non-zero provided 

lini( S 4) M', lim( S 4) T' and lim.( S 4) E' are finite and at least one is non-zero 
On both sides of the interface, the density reduces to 



(25) 







47T/i 



3 2 _ ■AMT't _ (6r) s / 3 M^ 3 E' 
2 M' 160M' 



if M' -> 0, 
otherwise, 



(26) 



and on both sides, 



X 
X 

R 
R 



/ - N 1/3 

, / rpl , M't , / 243T \ 171/ 



( 3M't 2 
\ 2M 



3T't 



E' 
10M 2 / 3 



I¥) 4/ T 



2r 
3r' 



(27) 
(28) 



Therefore \i and ^ and the shear are regular up to the junction. 



C. (c) Matching s3 to s5 



0. All the results in 



This junction is similar to (b). On approaching the junction in (s5), r\ ~ {<OTM) 1 '' i \E\ x l' 2 
section 5.2 for the kinematics and metric components follow with (s4) replaced by (s5). 

In this case R must increase in the opposite sense to that in section 5.2, i.e., here R must increase in the direction 
(s5) to (s3), as we will now show. Working in (s5), at the junction E = 0, and since E < 0, it follows that E' < 
in a neighbourhood of the junction. Also, near the junction the conditions (|23| ) hold, so R' and M' have the same 
sign and since M > and E may be as small as we please, R' and E' have the same sign. Since E' < 0, we have 
R' < 0, i.e., R increases in the direction from (s5) to (s3). It is interesting to note that, since R/R — > |er _1 > at 
a junction between (s3) and (s5), the continuity of R/R forces the existence of a finite region in (s5), adjoining the 
junction, where the azimuthal expansion rate R/R is positive even though all points in (s5) eventually satisfy R < 0. 



D. (d) Matching s4 to s5 

Both sides are constrained by E — > with M > and the result is obtained by combining results from (b) and (c) . 
In this case R must increase in the direction (s5) to (s4) by a similar argument to that given in (b). 

E. (e) Matching s5 to s6 

For solutions (s6), R' — 0, and so they may only be matched, across a comoving surface, to solutions of the type 
(s5), because R' — 0, X ^ requires E = — 1. This motivates a characterisation of (s6) within the family of LTB 
solutions by the conditions M' = T" = 0, M > and E = — 1. 

At this junction the (s6) side is unconstrained by the matching. The observer area distance R may increase in 
either direction on approaching the junction from (s5). Both metrics are regular on approach to the junction and it 
is less restrictive than the other four. 



When the conditions for matching and for a centre are combined, four classes emerge. 



VI. MODELS BY CLASS 

A. Open models with one centre 

By noting the sense in which R must increase at the interfaces (a) to (e) above, the only possible composite models 
are: 

C(sl) + , C(s2) + (s4) + , 

0(s2) + , 0(s3) + (s4) + , 

0(s3) + , C(s5) + 5(s5) + (s3) + , 

0(s4) + , 0(s5) + S(s5)+(s4) + , 

0(s5) + S, C(s5) + 5(s5) + (s3) + (s4) + , 

where O denotes a centre, and a superscript + (— ) implies that R increases (decreases) from left to right. Here S is 
any combination of (s5) , (s5) and (s6). Note that open models can be constructed from collapsing solutions [e.g. 
0(s5) (s6)]. Papapetrou |2(J] discussed a particular example of C(s5) (s3). 

In the above construction, we have noted from (0) that on t =const, d\ — \dR\j\/\ + E, where x ' IS radial proper 
distance. Hence by (pl|), if E > a > —1 for all \ > ft ( a i constants) then: 

JD 

X —* oo forces R — > oo if — — > 0, (29) 

dx 

d /? 

there is a finite value of y > 3 for which R = 0, if —r- < 0. (30) 

dx 

However, if _E ^ — 1 as^^oo, then neither of (|29|),(p0[) are necessary. 
An example of C(s5) (s5) _ in the class of open models with one centre is 

f _sr^j_ r 1 „ e -2r l f or <r<r , , Q1 , 

E = < sm 2 r L J u ' (31) 

\ -l + c- 2r for r>r 0l 

M = i S^ t M - + e_r °] for ° < r < ro ' (32) 

{Mcc+e-r for r>r , V ' 

T = 0, Moo > 0, 7T < r < 2tt, (33) 

and 

E= f -%[!--**] forO <r<r , 
I — 1 + c for r > tq , 

M= (^[M.-e--]forO<r<ro, 

[ Moo -e r for r > r , 

T = 0, < Moo < 2/3, r > 0, (36) 



is an example of C(s5) . In each o f (|33[) and (p6[), i? — >const> as x ~~ * °°- There are no spherically symmetric dust 
models with R — > as x — » oo [by (|2S|) and since, for the exact solutions, i? — > requires E 1 — > 0]. 

B. Open models with no centre 

By (p0[), to avoid a zero in i?, a model with no centre must either be composed entirely of (s6), or it must contain a 
section of (s5), in order to allow (at least one) minimum in R. Then the possible matchings are evident: 



(s3)-(s5)- II (s5)+(s3)+ 

(s4)-(s5)- r ) (s5)+(s4)+ 

( 8 4)-(s3)-(s5)- J { (s5) + (s3) + (s4)+ 



Examples and a detailed analysis of such models are provided in ||. In these models, due to the presence of collapsing 
solutions (s5),(s6), a centre does eventually form, but gravitational collapse will violate the regularity conditions in 
any case. 

C. Closed models with two centres 

These models must contain a region of (s5) , since there must be (at least one) turning point in R. The models cannot 
contain a region of (s2), (s4) or (s2)(s4), since the region would either contain a centre and match to another solution, 
or would match to other solutions on both sides. Hence E would vanish on both sides, and since E > throughout 
the domains of (s2) and (s4), E' could not have the same sign throughout, contrary to (p2) [with (Bj)]. There can 
be no (si) region in the closed model, since it does not match to any other solution. There can be no (s3) region in 
the model either, since R must increase in the direction (s5)— >(s3). Hence if (s3) contains a centre, it cannot match 
to (s5). Conversely, if (s3) does not contain a centre, it cannot match to (s5) on both sides, leaving the model open. 
This leaves just (s5) and (s6) to construct these models, and the possibilities are: 



D. Closed models with no centre 

Consider an SS dust model which has R > in some range < r < d (and at some t). This final possibility of 
composite models is obtained by identifying (matching) the surfaces r = and r = d. Since AR = 0, the model must 
be everywhere (s6) or else it must contain a region of (s5) [otherwise sign(i?') is constant in < r < d, which forces 
R(0) 7^ R(d)]. No regions composed from the solutions (sl)-(s4) may be present, since they would be forced to match 
to (s5) on both sides. This would force R' to change sign in the region (since R must increase away from (s5) into 



these solutions) and this is not possible, by (21). Hence the models may only be constructed from (s5) and (s6), with 
the possibilities: 

1S2 

where X denotes the surfaces which are identified (at which the standard matching conditions must be satisfied, as 
we have described). The spatial sections of these models have the topology of a 3-torus. An example is provided by 

ar 2 - 1 for < r < ±d, 

E = { a(r- \d) 2 - 1 for \d < r < \d> (37) 

a(r - d) 2 - 1 for |d < r < d, 

b + cr 2 for < r < \d, 

M=\ b+ \cd 2 - c(r - \d) 2 for \d < r < |d, (38) 

b + c(r - d) 2 for |d < r < d, 

1 J2 \ 4 



T = 0, a I 2b+-cd z I < -c, (39) 

where a, .., d are positive constants. Note that a closed model with no centre cannot be constructed from the homo- 
geneous (Friedmann-Lemaitre-Robertson- Walker) subclass of LTB (since the elliptic homogeneous solution has only 
one point with E = — 1, at which R is maximum). 

There are no further possible classes or composite models. There can be a number of different topologies constructed 
from these models but that involves different questions from those tackled here. Examples of models of types described 
in sections 6.1 to 6.3 are given in previous literature (see especially |9[]). 



VII. A BLACK HOLE IN AN EXPANDING UNIVERSE 

In this section we present an example which serves to illustrate some points of significance in using matched 
LTB models in cosmology. First that combinations of different LTB models can provide realistic exact models of 
cosmologically interesting phenomena; second that the richness of the models is often overlooked; thirdly that even 
in simple cases, there are subtle issues that need to be watched (for instance an apparently useful model may be 
flawed) and fourthly the illustration itself has some intrinsic interest as a model of the formation of a black hole in 
an Einstcin-de Sitter space-time which is different from the well known Oppenheimer-Snyder case. 

In his treatise on inhomogeneous cosmological models || Krasifiski gives the history of attempts to describe the 
formation of black holes in an expanding universe. The first successful attempt was by Barnes BJJ] , who proved suffi- 
cient conditions for a black hole to form and he studied several examples of collapse in LTB coordinates. Subsequent 
work on similar lines was done by Demiahski and Lasota |22|] and Polnarev p3 |. In 1978, Papapetrou published a 
model of a collapsing region in an Einstein-de Sitter spacetime. In his model the formation of the apparent horizon 
and the central singularity can easily be followed. In this model and in all others known to us, the singularity goes 
on accreting matter for an infinite LTB time, which limits their value in exact cosmology. In fact, although it is not 
widely known, this must happen irrespective of the precise nature of the collapsing part, unless it contains a vacuum 
region. The result follows from the matching conditions. If an s5 region is matched to either an s3 or s4 one then, 
as we have seen, E must vanish at the boundary. It follows that when r\ = 2-7T, r must be infinite on the boundary of 
the elliptic region (s5). Therefore, if the matter inside the elliptic region extends to its boundary then collapse to a 
singularity will take an infinite time to complete. We overcome this difficulty by arranging for the matter to extend 
only to some value of r within the elliptic (s5) region. 

Although the end result in our model is an Einstein-Strauss vacuole the dynamics of the formation of the hole is 
different. In Oppenheimer-Snyder collapse, the mass gets trapped at the boundary first and at the centre last. In our 
case, the reverse happens and the horizon begins to form at the centre and spreads outwards, and the singularity is 
covered at all times. 

A. The Model 

To define the model we choose the following arbitrary functions: 

T{r) = Vr, (40) 

B(r) = /-/3(>) 2 (l-£) 2 0<r<r o , (41) 

I r < r, 

\ar 3 tq < r, 



f 
M(r) = { 2 ar \ / 3 n<r<r , _ (42) 



V 3 fe)' 



r < r. 



The four arbitrary constants in the model are limited as follows: 

1. a> 

2. 16 > /3 > 

3. n > |r 

The following properties of the model follow easily from well known results: 

a. From a formula of Barnes |2lj , 

R' _ fM' E'\ (M 1 3E'\ tR 

~R ~ \1T ~ ~EJ ~ \~M ~ ~2EJ ~R' 

and the above, plus the appropriate exact solutions for R, it follows that R' > 0, as required to avoid a singularity 
in the metric. 

b. From condition (1) and R' > 0, the density is positive. 
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c. From (2) it follows that (1 + E) > 0, which ensures that r remains spacelike. 

d. The junction conditions are satisfied for the (combined) metric, and it is non-singular except at r\ = or 2-7T. 

For r > ro, the model is an Einstein-de Sitter universe with density p — (Girt 2 )^ 1 . For < r < 7*0, it represents an 
elliptic region which first expands and then collapses. In the region n < r < tq , it is vacuum, and for < r < r\ , it 
contains dust. 

B. Collapse 

Here we will discuss the dynamics of the elliptic region in some detail. It starts with a big bang at t = and all 
shells expand until they reach their maximum surface area which occurs at r\ = n, i.e. at times given by 

t= rEW (43) 

which is a monotonically increasing function of r. This means that shells with larger values of r reach their maximum 
later than those with smaller r. The shell bounding the matter, r — n, reaches its maximum at time 

t= 2[/3V2(i_ii ) ]3- I > 

After reaching their maximum surface area, the shells collapse to a singularity at R = 0, which occurs when r\ = 2tt, 
i.e. at time 

tt2M , x 

t = rEW 21 (45) 

which depends on r. The point r = is exceptional in that R vanishes for all t. However the behaviour of the density 
shows that after the universal singularity at t — 0, the point is non-singular until r\ = 2-k. In LTB time, the singularity 
at r — begins to form at time t(0) = to, where, taking limits, 

^ = -p£- (46) 

Note that the to here is twice that found by Papapetrou [E0[. The collapse process continues with the shells labelled 
r becoming singular at time 

t = t (l-^) ■ (47) 

The last shell at r = n becomes singular at t\ given by 

h =to(l-y\ , (48) 

after which the collapse is complete. 

The mass in the singularity at any time to < t < t\ is given by 

3 

(49) 



t \ V3 



1 r, 

M = -ari 

2 ° 

Therefore at the beginning of the collapse, M(to) = and at the end, M(t\) = ^arf. Given the size of a spherically 
symmetric tophat region and the mass enclosed, we could determine a. 

After time t\ the solution is only defined for r > r 1: i.e., in the exterior vacuum and Einstein-de Sitter regions. The 
vacuum represents the Schwarzschild region in comoving coordinates, [j 

For t > ti, the solution is equivalent to an Einstein-Strauss vacuole in an Einstein-de Sitter universe pi. 



2 



Singularities continue to form for r > n at times 



m 



3 / r \ -3 



^ ro / V ro 
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C. Horizons 

As usual in LTB spacetimes, we use the apparent horizon as diagnostic for the existence of an event horizon and 
therefore a black hole. The apparent horizon is given by [2l]]P| 

R(r,t)=2M(r), R<0. (50) 

This formula can be put in the alternative forms 

R=-{l = E) 1 ' 2 (51) 

and 

sin(ry/2) = (-£T) 1/2 , tt < r, < 2tt. (52) 

To simplify the notation, we define w := ~. Then ( |5^ ) gives 

sin(ry/2) = (3 1,2 w(l -w), n < rj < 2tt, (53) 

and from the definition of r) in solution (s5) and the definitions of the arbitrary functions, 

toff) — sin 77) 

tAH(r) = -^j- #, 0<w<w u 54 

27r(l — wy 



tAH 



/ n (t (r) - smrj)\ /wi\ 



where w\ := ri/ro and 77 is determined by (|53|). Equations (54) and (|55| ) together give the equation of the apparent 
horizon. The areal radius of the apparent horizon at coordinate r is given by 

R AH = ar 3 , 0<r<r x , (56) 

Rah = a.r\, r\ < r < r , (57) 

where we have used the definition of M for the appropriate range of r and equation ( p0[ ) . 

From (B3h if w = then 7/ = 27r and so (p3) gives tAH — to- Thus at r = the singularity and the apparent horizon 
form together. When this happens, the singularity may be naked Pq-p8||, but we will not discuss that here. We will 
concentrate instead on the formation of the black hole and assume that < w < 1. 

We denote the time at which the singularity forms at r by t s (r), which is given by (|47|), i.e., 

t.(r) = to(l - w)- 3 . (58) 

Then for < r < n, i.e., inside the collapsing dust sphere , we obtain from (|5 

io(27r — r\ + sin 77) 
2tt(1 - w) 
and for r\ < r < ro, i.e., in the surrounding vacuum, 



t s (r) - tAH(H = — ^37^ — ^3 , (59) 



. / x , /s t (27T-7/ + sin7/) ( Wl \?> 

t s {r)-t AH {r)= Ml _ w)3 (-) , (60) 

where 77 is given by equation (p3). Given the range of 77, t s — t A H > always. So the apparent horizon forms first 
and the singularity is not naked. This agrees with a result of Joshi |29) . The horizons starts to form at the centre 
and spreads outwards with time. This is different from the behaviour of the Oppenheimer-Snyder solution for which 
the boundary and the mass get trapped first and the centre last. 

In the limiting case where w — ► 1, the inner collapsing sphere extends to the Einstein-de Sitter exterior, so that 
there is no vacuum region. This is the case considered by Papapetrou and analysis of the limits confirms his result. 
Both t s and t A H tend to 00. However the difference between them remains finite and t s remains greater than t A n- 



and are to be interpreted as the arrival of successive shells of test particles which label the coordinates of the vacuum in this 
gauge. This process continues until the shell r — ro arrives at time t — 00. 

3 This is the apparent horizon bounding trapped surfaces associated with the collapse, which begins at 77 = ir. There is also 
an apparent horizon related to the initial expansion which all particles must cross during their expansion phase from the big 
bang white hole. Gautreau and Cohen [E5| call this a boundary of expelled surfaces. 
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VIII. CONCLUSION 

Here a set of spherically symmetric inhomogeneous dust models has been provided which can be used to construct 
cosmological solutions to Einstein's field equations for a range of astrophysical situations including voids and black 
holes. A section of the book by Krazihski H is devoted to applications of these models in cosmology. In setting 
out the models new results have been obtained which fill gaps in the literature. In particular we have demonstrated 
the matching of the Kantowski-Sachs solutions and derived a new solution where the spatial sections have a torus 
topology. Also the particular regularity conditions we use lead to a restricted but, we would argue, cosmologically 
more useful set of matched solutions. The categorisation we provide of the allowed cases with centres is useful for the 
construction of models. 

The value of these matched solutions is demonstrated in the final example of black hole formation which illustrates 
very clearly why it is important to consider details when using inhomogeneous models. We have shown that using 
matched LTB solutions a model can be constructed to describe collapse to a black hole in an Einstein de Sitter 
background but it must contain a vacuum region if the singularity is not to continue to accrete matter for an infinite 
LTB time. 
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